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in 

We present the ground state wave functions for systems of one- dimensional interacting fermions. 
It is shown that these systems undergo phase transitions similar to the Kosterlitz-Thouless one 
\ independently of the interaction details. In the limit of an infinitely strong interaction the phase 

transition turns into the usual second order phase transition in a chiral phase. The temperature of 
■ the phase transition is calculated. 
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I. INTRODUCTION 



One dimensional fermion systems were a subject of intensive studies both in field theory and in condensed matter 
physics. As was demonstrated by Tomonaga Q and Luttinger in their pioneering papers the long range excitations 
of such a system (under rather general conditions) can be expressed in terms of non-interacting bosons. These degrees 
of freedom were made explicit in the elegant method of bosonization proposed by Mattis and Lecb 3] . The recent 
interest in this field is mainly due to the development of submicron techniques which allowed to produce very pure 
quantum wires. In such wires only few levels (or sometimes even one) corresponding to the quantization of electrons 



, in perpendicular directions are occupied. In spite of the fact that the number of the experimental papers in the field 
is still not numerous (see, e.g. 3, 0], 0) the systems under discussion are, in principle, accessible by experiment. 
The bosonization technique allows one to calculate all correlation functions of the "density-density" type for systems 
i ^ i . of interacting fermions in one dimension. However it says nothing about the ground state of this system. The 
correlation functions reveal a number of anomalies of the fermion system (see @, 0]j E3) : they have oscillating 
contributions with wave vectors equal to 2pf or Apf which decay very slowly with distance. In the literature these 
contributions were interpreted as follows: the oscillations with Fermi momentum, pf, doubled were related to the 
Peierls instability (connected with the charge density wave [Tll |. 0) while the oscillations with 4p/ frequency were 
interpreted as a marginal Wigner crystal |12(. It is commonly believed that the system under discussion is a kind 
of normal liquid, because quantum fluctuations destroy any order parameter, and a phase with long-range order is 
impossible even in the zero temperature region The common point of view was put into the words: "Luttinger 
' liquid is a normal (not symmetry-broken) metallic phase" |13| with gapless boson spectrum. This viewpoint is us ually 
supported by the Landau theorem |14| which states that in one dimension a long-range order cannot exist |15| . 
' However two more points should be taken into account: 

a ; 

• It is well-known that the Landau theorem is app licable not to all systems in low dimensions. Some of them 
I ■ appear to be in the Kosterlitz-Thouless phase [lfj . In such a phase the order parameter density tends to zero in 

an infinite system. However a long-range order exists because correlation functions decay as some power of the 
distance. This means that the correlation is present in a whole specimen. As a result the macroscopic properties 
of the system are quite similar to the system with broken symmetry. A Kosterlitz-Thouless phase can form if 
the gapless excitations, which should be present in the system after spontaneous symmetry breaking due to the 
Goldstone theorem ([lfl, 0), do not interact. 

At non-zero temperatures (0 7^ 0) this phase can appear in two-dimensional systems; at Q — a Kosterlitz- 
Thouless phase is possible in one dimension. (The crossover temperature and its relation to the phase transition 
temperature <d c will be considered below, see section lllT|) . 

• The Goldstone theorem itself can be broken in one dimensional models with Adler-Schwinger anomaly. For that 
to happen the interaction has to be strong enough. As an illustration one can consider the massless Schwinger 
model [l^. I n the Coulomb gauge this model is a particular case of the Luttinger model where the potential 
is a linear function of the distance between electrons. This gives rise to the Goldstone theorem violation: all 
excitations have a gap in spite of spontaneous breaking of the chiral symmetry. For this reason, excitations 
cannot suppress the long-range order at temperatures smaller than the gap. As a result one has a second order 
phase transition in one dimension even at finite temperature |20j. 

This makes the statement that the Luttinger liquid at = is in an unbroken phase doubtful. 
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In order to clarify this question we calculate the wave functions of the ground states in the Tomonaga-Luttinger 
model. While the details of these wave functions could depend on the interacting potential, qualitatively all possible 
ground states reveal the same phenomenon. 

In one dimension the Fermi surface reduces to two isolated points in phase space (p = ±p/)- One can neglect 
transitions between these two points. This is a good approximation, at least if the potential is a decreasing function 
of the momentum transfer. As a result, the number of electrons near every point (left- and rig/it-particles) should 
be conserved and the system acquires complementary (chiral) symmetry. It is this symmetry which, as we will see, 
breaks down spontaneously in the model. 

The Fermi surface in one dimension is absolutely unstable: even infinitesimal electron-electron interaction leads 
to instability of the electron ground state. The electron distribution function is of the order of 1/2 near the Fermi 
surface [2]]]. This means that there is a hole located near every electron. Naturally, they attract each other and form 
a kind of bound state consisting of a right electron and a left hole (i?Z-pair). This is quite similar to the formation 
of a Cooper pair in a superconductor but the quantum numbers of the bound state are different: instead of non-zero 
electric charge it has non-zero chirality p^ . 

Of course, this fact itself is not enough to speak about a new correlated phase. We will check using the explicit 
ground state wave function constructed in the paper that long-range order is indeed present in the system. At low 
temperatures Luttinger liquid undergoes a phase transition similar to the Kosterlitz-Thouless one. 

We will see that in the limit of the infinitely strong interaction and at zero temperature the system is in a phase 
with broken chiral symmetry and non-vanishing order parameter. The properties of the Luttinger liquid in this limit 
are quite analogous to the properties of the massless Schwinger model where the spontaneous breakdown of the chiral 
symmetry is well-known. This can be expected as the interaction in the Schwinger model is also infinite (it grows 
with distance). On the other hand, in contrast to the Schwinger model, the spectrum of the Luttinger model remains 
gapless. The only effect of the interaction on the spectrum is the renormalization of the Fermi velocity v'j which is 
proportional to the strength of the interaction. 

If the interaction in the Luttinger model is considered to be finite, the order parameter vanishes (in the infinite 
system). The system appears to be in a Kosterlitz-Thouless phase with correlators which decay as some power of the 
distance. However the properties of this phase are rather close to those for a system with non-zero order parameter. 
The Kosterlitz-Thouless system transforms smoothly to the phase with non-zero order parameter as the coupling 
constant increases. 

Another matter is the temperature O c of the phase transition to the unbroken phase. This temperature is always 
of the order of the gap in the spectrum. For this reason O c is finite for the Schwinger model and of order of v'j/ L 
for the Luttinger liquid (where L is the size of the system). In other words O c tends to zero for an infinite system. 
However for real systems O c is not so small. Indeed, if one takes Vf ~ 10 7 cm/ sec and L ~ 10 _4 cm then 9 C ~ 1° K 
■Vf/vf. This means that the theory based on the limit L — > 00 can be inadequate. Of course, in order to discuss O c 
for finite size systems we have to define the notion of phase transition in this case. This is done in Section IlIII 

We will see that a non-zero order parameter is the reason of the anomalies in the correlation functions mentioned 
above. Anomalies exist even at temperatures higher than O c owing to fluctuations of the chiral phase in the phase 
with non-broken chiral symmetry. There is a clear-cut distinction between the chiral phase and the Peierls's phase. 
The latter is the second order phase transition in the phonon system and the chiral symmetry of electrons is broken 
explicitly, in the Hamiltonian. On the contrary the chiral transition originates from the spontaneous symmetry 
breaking in the electron system. 

Finally, let us discuss effects which may indicate the existence of a condensate. Obviously a charge-neutral con- 
densate can not reveal itself in experiments associated with charge transfer. However it contributes to the effects 
concerned with energy currents and should not transfer heat. Hence one can think about thermal anomalies related 
to the condensate. In fact, we keep in mind the effect similar to the thermomechanical effect in supcrfluid helium. 
(The temperature decreases with increase of the superfluid mass |2^|). We plan to discuss this problem in a separate 
publication. 

The paper is organized as follows. In Section^we introduce the Hamiltonian, definitions of left- and right-particles 
and so on. We present our results and discussion in Section ITTT1 We relegate the derivation to the next sections as 
the calculation is a bit cumbersome. 
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II. NOTATIONS AND GENERAL EQUATIONS. 

We begin with the usual expression for the Hamiltonian of interacting spinless electrons in one dimension: 

H = j dx&(x) (-±.%-^%(x) 
dxdy& (x) # (a;) V (x - y) & (y) * (y) . 



(1) 



Here $ (a;) is the electron field, m is the electron mass and p, is the chemical potential. V (x — y) is the electron-electron 
interaction which we will discuss below (K = 1). 

As usual let us separate in the electron wave functions left- and right- particles p4|: 



^ (x) — exp (ipfx) i&R (x) + exp (—ipfx) $l (i). 



(2) 



It is implied here that the wave functions wjj £ (x) are varying over distances much larger than l/pf. We also restrict 
ourselves to the Tomonaga-Luttinger model and assume that the characteristic scale of the potential is large as 
compared to 1/p/. Lastly, for simplicity we consider only an electrically neutral system where the a positive charge 
of ions is distributed homogeneously along the channel. 

Let us proceed to the electron-hole representation for the right (left) particles: 



* ( •'• ) = / y~ ( exp ( ±j -P x ) ®R,L (?) + ex P (Tipx) b RL (pfj = d R . L (x) + b RL (x) . (3) 



Here fit (a) and b'(b) are creation (annihilation) operators for electrons and holes. The Hamiltonian of eq.Q can be 
written for neutral system in terms of electron density operator 



in the form 



H = dx 



g(x) = g R (x) + ql (x) 



y R (x) v f (-id x ) (x) + *t (x) v f id x y L (x) 



dxdyg (x) V (x - y) g (y) . 



Here Vf is the Fermi velocity, and 



Qr,l (x) = d R L (x) a R x (x) - b R L (x) b R>L (x) + d R L (x) b R L (x) + b RnL (x) a R , L (x) . 



(4) 



(5) 



The expression 0| coincides with the Hamiltonian of the Tomonaga-Luttinger model (see, e.g., Ref. |24j). 

Let us discuss now the electron-electron interaction V(x — y). Its form depends on the relation between the usual 
3D screening radius Rd (for simplicity we will consider the case of Debye screening) and the transverse size of the 
channel d. Indeed one has to take into account that the electrons are one-dimensional only for distances \x — y\ which 
are much larger than d. Therefore if Rd <C d one can use a point-like interaction Vb<5 (x — v) |25| . This is the case 
for metals |2q . In the opposite case Rd ~S> d, semiconductor, one should use an ordinary Coulomb potential. In one 
dimension the Fourier transform of the Coulomb potential is logarithmically divergent. This divergence is regularized 
by the smaller of two quantities, either pf or the inverse channel size, Thus |27| 



• At R D < d 



• At R D > d 



V(p) = V , 



2e 2 logf|f) forl/ P/ »d 
2e2log (w0 for 1/p/ «d 



V(p) = 



(6) 



(7) 
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The Hamiltonian of the Luttinger model, eq.QJ, presented in terms of electrons and holes is completely defined 
without any additional regularization of electron operators. In particular, the commutator of the R- and L-densities 
in this representation reproduces the well-known Schwinger anomaly [lflj : 

i d 

[qr,l (x) , Qr,l (y)\ = ± '^g^ S i x ~y)- ( 8 ) 

These relations are the starting point of the bosonization technique. Usually one derives eq.© regularizing the 
product of ^-operators by a small shift of their arguments |2^. This is not necessary, however, in the electron-hole 
representation |2(j as the creation-annihilation operators for R,L electrons (a^ L (x)) and holes (b^ L (x)) are non-local 
in the coordinate space: 

{a R (x) ,a R (xi)} = {b^ (x) , b R (xi)} = ^- ■ 1 _ „ (9) 



{a{ (x) , a L ( Xl )} = {b{ (x) , b L (a*)} = — ■ Xi _ x _ iS - (10) 

(in momentum space these anticommutators are (5-functions). Using these anticommutators for the densities of right 
and left electrons Qr^l (x) in the form of eq.10 we immediately reproduce the Schwinger anomaly. This means 
that, being formulated in the electron-hole representation, our theory is completely determined without any further 
redefinition of the density operators. 

The Hamiltonian 0] is invariant both under vector transformations: 

= e iav * R , * L (x) = e iav * L (11) 

and under chiral transformations 

9 R (x) = e ia °* a , * L (x)=e- ia °* L . (12) 

The first invariance leads to conservation of electrical charge (the number of left electrons plus the number of right 
ones), the second means that chiral charge (the difference of numbers of right and left electrons) is conserved. However 
we will see below that the ground state of the model is constructed in such a way that the second symmetry can be 
spontaneously broken. 



III. APPROACH, RESULTS AND DISCUSSIONS 

The standard approach to systems of many particles is based on Green functions. The one-particle Green function 
gives the information about the spectrum of excitations; the many-particle Green functions allow one to calculate 
different correlation and response functions. Of course, the Green functions give some information about the wave 
functions of the states but this information is indirect. 

In principle, one can obtain the wave functions of stationary states (and, in particular, of the ground state) by 
solving the corresponding Schrodinger equation directly. However for systems with infinite number of degrees of 
freedom this equation is too complicated. A more practical approach can be based on the evolution operator [2jJ 
which is a sum : 

S (T) = |n >< n\ exp -iHT\m >< m\. (13) 

m,n 

Here \n > are the exact wave functions of the Hamiltonian H in the second quantized representation, T is the time 
of observation. The evolution operator determines the evolution of an arbitrary initial wave function (< m|) from the 
time t = up to final states |n > (at t — T). (We imply from now that the Schrodinger representation for operators 
with time-dependent wave-functions is used.) 

Formula suggests the general method to obtain wave functions. One has to calculate first the evolution operator 
and represent it as a sum of time-dependent exponents. The coefficients in front of these exponentials are products of 
exact wave functions and their complex conjugates. In order to extract the ground state wave function one has to take 
the limit T — > oo (we add an infinitesimal imaginary part to the energy). Proceeding to Euclidean time (T — > —i/Q 
we see that evolution operator determines the density matrix for the equilibrium system at non-zero temperature (see 
end of the Section TlV All . 
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The advantage of this method is that the evolution operator can be written explicitly as a functional integral with 
definite boundary conditions (see ea. l|26[l )• This functional integral is rather simple for the Luttinger model with 
Hamiltonian (@J and can be calculated exactly. This allows one to construct wave functions of all states in the model 
and, in particular, the ground state. This will be quite enough in order to demonstrate the symmetry breaking and 
to calculate the temperature of the phase transition. 

We will keep the size of the system finite. This is important not only to regularize infrared divergences in the system 
but mostly because of the peculiar situation with the temperature O c mentioned in the Introduction. However, first, 
we have to define the concept of phase transition in finite size systems. 

Usually, the critical temperature is defined as a point where thermodynamic quantities have a singularity. Of course, 
this is the case only in a infinite system because all singularities smear out if the size of the system is finite. The same 
is true for the coherence length — it cannot be larger than the size of the system. 

In this paper we will adopt the point of view suggested by Landau in order to describe the second order phase 
transitions [30(. He introduced the order parameter as the main quantity for the description of phase transitions 
related to the spontaneous symmetry breakdown. By definition, the order parameter is zero in the high symmetry 
phase (with the same symmetry as the Hamiltonian) and non-zero in the phase with broken symmetry. In the case 
of chiral symmetry ea. l|12|) the following quantity can serve as an order parameter: 



This quantity is not invariant under transformations of ea. (|12fl and should be zero if the chiral symmetry remains 
unbroken. Note that we use a macroscopic order parameter (integral over the whole size of the system). In the broken 
phase this quantity is proportional to the volume of the system. 

The Kosterlitz-Thouless phase represents the intermediate case when 



so A is still infinite in the thermodynamic limit while the density or order parameter A/L vanishes. Let us point 
out that A appears to be non-zero even at O > G c due to fluctuations of the broken phase in the phase with higher 
symmetry. What is important for the latter case is the fact that A does not increase with L. 

Intensive thermodynamic quantities remain smooth for finite size system even at the point of the phase transition. 
However what matters is the fact that they depend explicitly on the system size and tend to infinity (or acquire a 
jump) at L — ► oo. 

Usually one proves that the system is in the Kosterlitz-Thouless phase by investigating the behaviour of the four- 
fermion correlator which does not break the chiral invariance ( below in section flV A lt ve will consider such a correlator, 
namely the probability to find an i?L-pair at a large distance r from the LR pair). If such a correlator decreases 
sufficiently slowly with the distance, the system is in the Kosterlitz-Thouless phase. The limiting case when the 
correlator remains constant at large distances corresponds to a non-zero density of the order parameter and ordinary 
broken symmetry. In fact, this definition of the Kosterlitz-Thouless phase is equivalent to our definition given above 
but definition ^141 115(1 is more convenient for us. 

In one dimension the Kosterlitz-Thouless phase can exist only at O = or, to be more precise, for the temperatures 
which tends to zero at L — > oo. There is no need in microscopic theory in order to estimate C : one can use general 
phenomenology applicable to all Kosterlitz-Thouless systems (see, e.g. |2^). 

Let us assume that the chiral symmetry is indeed spontaneously broken in the Luttinger model (of course, this can 
be proved only in microscopic theory). According to the Goldstone theorem the chiral phase a c (the phase of the 
^Jj^l operator) becomes a massless boson field. In the long-range limit only fluctuations of this field are relevant 
and its effective Lagrangian reduces to (Euclidean time r = it is used, as we want to consider non-zero temperatures 
below): 



where V and w are phenomenological constants (calculable in microscopic theory). 

To judge if the system is in the Kosterlitz-Thouless phase it is enough to consider the behaviour at large distances 
\x — y\ of the chirality conserving correlator: 




(14) 



< a T < 1, 



(15) 




(16) 





(18) 
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Calculating the latter integral at = we get: 



F(x — y) ~ exp 



2V 



dk dk sin 2 [l/2(k,x-y) 



(27T) 



7 2 k 2 ) 



(19) 



The two-dimensional integral ( one space one time dimension ) in the exponent Ijl9(l diverges logarithmically and 
hence 



F(x - y) 



V 2 



{k m ax\ x y\) 



1/2ttV 2 w ' 



(20) 



This proves the existence of Kosterlitz-Thouless phase at O = 0. 

If the temperature is non-zero the integral over kg in ea. (|19|l should be replaced by a sum over discrete values 
kg = 27m0 (n integer). At high temperatures only the term with n = survives at large distances and we are left 
with an one-dimensional integral in k which leads to the correlator exponentially decreasing with distance: 



F(x — y) ~ V 2 exp 



e 



2kV 2 w 2 



\x - y\ 



(21) 



Clearly this correlator describes the unbroken phase. 

A power - like behavior of the correlator of l|20|) is valid in the region 



For 



\x — y\ < w/Q. 



e < e c = w/l 



this takes place for the whole specimen, i.e. the system is in broken a phase. The temperature O c is a temperature 
of the phase transition. 

One can recognize in this estimate the excitation energy with the smallest momentum possible in a finite size system. 
In the Luttinger model this energy is equal to uj m in — 27rwJ / L with renormalized Fermi velocity vj — Vf y/ 1 + Vq /irvf 
0- We see that if the spectrum of excitations is gapless (like in the Luttinger model), then the temperature of phase 
transition is inversely proportional to the specimen length. This result can be obtained in the microscopic theory as 
well (see Section TlV All . 

Turning to the microscopical theory we begin with the simplest case: the short range potential in the limit of an 
infinitely strong interaction © 



TTVf_ 



< 1. 



(22) 



In the leading order in this parameter the evolution operator appears to be very simple and the wave function of the 
ground state can be represented in close form. In the temperature region 



Q chiral 

the ground state wave function is of the form 



Zq exp 



dx exp {i0) at, (a;) b L (x) 



dyexp{-id) at (y) bt (y) 



\F>. 



(23) 



(24) 



Here \F > is the filled Fermi sphere and Zq is the normalization coefficient. There is an infinite set of degenerate 
ground states labeled by the continuous parameter 9 which has the meaning of the order parameter phase. 

One can see the symmetry breaking immediately because the wave function (|24|l is not invariant under the chiral 
transformation l|12|l . Besides one can check directly that A oc L [22|. It means that there is a second-order phase 
transition in this limit. 

The wave function of l|24l) is a mixture of states with different chirality. (We assign chirality +1 to a right electron 
and a left hole and —1 to their counterparts. So bosons in ea. l|24|l are neutral in terms of electric charge but have 
a non-zero chirality, ±2 ). Such a ground state implies that the states with different chirality are all degenerate in 
energy. This is a typical situation for systems with condensate: the addition of one pair to the condensate does not 
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cost any energy. However this degeneracy is possible only if the size of the system is large enough, namely, we will see 
that it should be L » L min ~ 2irvf/<d. At L < L min the ground state has fixed chirality (equal to zero) (see ea.i|62|)) 
and the order parameter A vanishes, i.e. there is no spontaneous symmetry breaking. These considerations put a 
lower bound on the temperature region where a chiral phase can exist: ^> ® chiral- So, ® chiral is the degeneration 
temperature. 

Let us estimate the density of chiral pairs in the ground state. The wave function (|24|) implies that all electrons 
are bound to the pairs. Hence the density of RL coincides with the density of i?-electrons: 

N R (p) = < n|4 (p) a R (p) \Q > e = L/2. (25) 

(see H2)- This quantity reflects a well-known fact: the distribution function of electrons is of the order of 1/2 near 
the Fermi surface |2l|. If the interaction is infinitely large all electrons and all holes are bound in exciton-like pairs. 
As a result we get the value 1(251) which is maximal possible. 

In the model under consideration Nr (p) is momentum independent and the total number of pairs Nr diverges at 
large p. (This is the defect of the point-like electron-electron interaction.) The sum over all states should be restricted 
either by pf or, at pfd 3> 1, by the inverse size of the channel because at larger p electrons cannot be considered as 
one - dimensional (see |26J). As a result 




Thus the number of pairs Nr is only a small fraction of the total number of electrons (Lpf/2n). This does not 
mean, of course, that in this case the Luttinger liquid behaves like a normal one. The reaction of the system to the 
slowly varying external fields is determined completely by the electrons near the Fermi surface which are all paired. 
This situation reminds superfluid helium where (even at zero temperature) the density of the condensate is only few 
percents of the total one. Nevertheless the whole mass of helium is superfluid |33| . 

Let us proceed with the region of high temperatures: 6>8 C = ^>min- I n this region the macroscopic order param- 
eter (|14f) is proportional not to the volume of the system but to the some characteristic length £ (0) = v c J (0 — C ) 
(see the end of section llV A|) and the density of the order parameter A/L vanishes in the limit L — > oo as it should 
be. Hence the temperature C has indeed the meaning of the temperature of phase transition from the symmetrical 
phase to the phase with broken chiral symmetry. 

The length £ (0) plays the role of coherency length in our system. At lengths less than £ (0) the wave function of 
the system coincides with the coherent exponent (12411 . However at larger distances the order disappears. 

The macroscopic order parameter A can be non-zero even in the symmetrical phase due to fluctuations of the 
broken to the unbroken phase. What matters is the behaviour of A with the size of the system. If A does not increase 
with L (A oc £ with £ finite), we deal with unbroken phase, where A increases with L the long range order appears 
(and £ ~ L). Such a dependence C on L can be considered as the definition of symmetry breakdown for a finite size 
system too. 

On the other hand it is obvious from such a definition that a temperature of the phase transition in the finite size 
systems can be defined only up to l/L corrections and the phase transition is smooth within the l/L region around 
the temperature of the phase transition. In the Luttinger model, where the temperature C itself is of l/L order, we 
can define C only up to a factor of order of unity. This is the price we have to pay for considering a phase transition 
of a large but finite size system. However there is still a clear-cut distinction between the case with a correlation 
length of the order of the size of the system (broken phase) and the case when £ <C L (unbroken phase). 

As was already pointed out in the Introduction, the case of the infinitely strong interaction is very special. We 
will see that if the interaction is finite, then the macroscopic order parameter A grows with the system size but more 
slowly than L (at <C C ). In the case of a short range potential (see Section liV A 1(1 A behaves as some power of L. 
This corresponds literally to the definition of the Kosterlitz-Thouless phase. If we consider the potential of Coulomb 
type, then A depends on L in a more complicated way fsection flV C|l but still A increases with the size L. Physically, 
this case is quite similar to the usual Kosterlitz-Thouless one. 

To summarize, the Luttinger model at < C is always in Kosterlitz-Thouless phase with broken chiral symmetry. 
At ~ C it undergoes a phase transition which in the limit of infinite interaction turns into the ordinary second 
order phase transition. 

IV. GROUND STATE OF THE TOMONAGA-LUTTINGER MODEL. 

The evolution operator 1(13(1 of the quantum system can be represented as a functional integral with definite boundary 
conditions (see, e.g. |2j|). Usually one derives this representation for boson systems, for the sake of completeness we 
give in Appendix lAl the derivation for fermions. 
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The theory with arbitrary electron-electron interaction can be reduced to a theory in an external field by means of 
the Hubbard-Stratonovich transformation (34[(see below, ea. (|37[l ). One has to integrate over the value of the external 
field in order to return to the original 4-fermion interaction. For this reason we consider first the evolution operator 
for one-dimensional electrons placed into an external field Q(x,t). It is of the following form: 

§($)=[_ V^V^cxp S (*,*). (26) 
Here fy, ^ is the electron field (Grassmann variables) and S is the action: 

S = i J dt J dxH R (x, t) [id t - v f id x + $ (x, t)] ^r (x, t) 

+ (R,v f ~L,-v f ,). (27) 

Integration over 'J, \1/ in ea. l|2tj|> is performed with given boundary conditions at t = and t = T: 
At t -> +0 

^r,l {x, t) = clr.l (x) + arbitrary negative frequency part 
^r,l {x, t) = b Rt L (x) + arbitrary negative frequency part 

At t -> T - 

^r,l (x, t) — b]^ L (x) + arbitrary positive frequency part 



^r,l (x, t) = a' R L (x) + arbitrary positive frequency part (28) 

The creation operators of electrons and holes d + , b + are the variables which enter the wave functions of the states 
in the sum of ea. (|13(l . The annihilation operators a,b enter the conjugate wave functions. They anticommutc: 
{a, a + } — {b , b + } = since they belong to different instances of time as long as one calculates the evolution operator 
(see Appendix lAl for details). 

It is possible to separate explicitly the dependence on creation-annihilation operators for the evolution operator in 
a given external field determined by the functional integral l|2t)|) . Let us introduce new integration variables: 

■*R,L = ^R.L + XR,L 

*R,l = *r,l + Xr,l- (29) 

The saddle-point fields l are supposed to obey the Schrodinger equation in the external field Q(x,t) with 
given boundary conditions l|28|) . The "quantum" fields XRL(x,t) are arbitrary but obey zero boundary conditions: 
\,;.,r/.: 0. 

The solutions ty R L can be represented in terms of the Feynman Green function in the finite time Gr^l which is 
defined as follows. It is a solution of the Schrodinger equation: 

[id t T v f id x + $ (x, t)} Grx (x, t; x x ,h) = ^ (2) (x-x u t- h) (30) 

with the following boundary conditions: at t — > +0 the Green function Gr(x, t, x±, ti) should coincide with the Green 
function of free fermions in the lower semiplane of the complex variable x (being arbitrary in the upper semiplane) . At 
t — ► (T — 0) it coincides with the free Green function in the upper semiplane. For the Green function of left electrons 
GL{x,t,xi,ti) one has to exchange upper and lower semiplanes. 
The free Feynman Green function is equal to @: 

G° R L (x,t;xx,tx) = [vf (t - h) T{x- xi) - i<5sign(i - ii)] -1 . (31) 

ZTTl 

In one dimension the Schrodinger equation (|30|) can be solved for an arbitrary external field &(x,t): 

Gr, l (x, t; X!,t\) = G° RL (x, t;x!,ti)x 
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x exp 



i dt' dy^{y,t')(G% L {x 1 t- V X)-G% L {x ll t 1 \y,t 1 )) 



(32) 



Now it is easy to verify that the saddle point fields ^> R L can be expressed in terms of these Green functions as 
follows: 



*fl l (x, t)= dx' \Gr, l (x, t; x', 0) a R , L (x') - Gr.l (x, t; x'T) b RL [x') 



*°r,l (x, t) = - dx' G R>L (x', 0; x, t) b RtL (x') - G R . L (x', T; x, t) a R L (x r ) 



(33) 



In order to check that these fields obey the required boundary conditions let us note that d R (x) and b R (x) are regular 
in the upper semiplane (see eq.lJSJ). Therefore the positive frequency part at t — > +0 of G R (x, t, x%, ti) is determined 
by the pole contribution at x 1 = x + iS and is equal to a R (x) as it should be. The second term in ea. (|33|l gives 
the negative frequency part which is arbitrary. Similarly, one can check the boundary condition also at t — > (T — 0). 
Inside the time interval (0, T) the saddle point fields satisfy the Schrodinger equation as is seen from ea . (|3()|l for the 
Green functions. 

The contribution of the saddle-point field to the action is: 



i=R,L 



Sq = I dxdx' bi(x')Gi(x' ,0;x,e)di(x) + d\ (x')Gi(x' ,T;x,T — e)b\ (x) 

- a\ (x') Gi (x', T; x, 0) a t (x) - k (x') G t (x', 0; x, T) b\ (x) 



(34) 



Since the saddle point fields obey the Schrodinger equation there is no term linear in the quantum field \ in the 
action. 

Dependence of the evolution operator in the external field on the creation-annihilation fcrmion operators is com- 
pletely determined by ea. H34|) . Integral over quantum fluctuations produces the determinant of the Schrodinger 
operator in the external field $ (it is calculated in the Appendix [5J| : 



1 f^ 1 f°° dv 

log[Det$(T)]=-—J dtdhj ^$(-p,t)$(p ) t 1 )|p|expHH«/l*-*i 

The complete expression for the evolution operator in the external field has the form: 

S ($) = exp (S + log[ Det<S> (T) })\F >< F\, 



(35) 



(36) 



Now we can express the evolution operator for the system of interacting fermions in terms of this operator. We will 
use the well-known identity |34j : 



exp 



i 

2 Jo 



dt 



I 

2 Jo 



i 

2 

dp 
2^ 



dt 



dt 



V(p)g{p,t)g(-p,t) 



dp 
2^ 

^{p^^i-p^V-'ip) 



(g(p,t)^(-p,t) + g(-p,t)^(p, t)) 



Here V (p) is the Fourier transform of the interacting potential. The normalization coefficient J\f is 



TV = / 2?$ exp 



J J — CXD 



(37) 



(38) 



In order to prove ea. (|37|) it is sufficient to shift the variable of integration $ to $ — V g in the integral 

Z?$ exp 



/ 



-Tr ($t$y-l) 
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Applying identity (|37[) to the functional integral that determines the evolution operator for the Tomonaga-Luttinger 
model we express it in terms of the evolution operator in the external field at the price of an additional functional 
integration over the scalar field t): 



S e - e = 1 J V<S> exp 1 £ dt J°° (p, i) $ (-p, i) F- 1 (p) 



(39) 



Expression (|39() is explicit: while it is not possible to perform the final integration in $(a;,t) in closed form, it is 
easy to obtain an arbitrary term of the evolution operator expanding it in creation-annihilation operators. This will 
be enough in order to calculate the evolution operator. 

Indeed, let us expand the evolution operator in the external field in powers Sq. The arbitrary term of the expansion 
contains a number of Green functions in the external field (|32|l which are exponents linear in the external field. Together 
with the action 1)39(1 and the determinant, ea . I(38|l . we get an integral of Gaussian type in &(x,t) which can be easily 
performed. The result of the integration depends of the electron-hole configuration considered. 

Let us introduce the following system of notations for the coordinates entering the electron-hole creation-annihilation 
operators: 

1. We will denote by x the coordinates of the right particles and by y the coordinates of the left ones. 

2. We will put the tilde on coordinates related to annihilation operators (initial state) and leave coordinates of 
creation operators (final state) without tilde. 

3. We will put primes on coordinates which are related to holes. 

It is convenient to proceed in the exponents of the Green functions (|32|) to momentum space using the expression for 
the free Feynman Green functions: 

Gr, l (p, Mi) = S± P 9 (t-h) exp [Tipv f (t -h)]- 9 TP 9 {h - t) exp feipv f (t -h)}. (40) 

Collecting all terms in the exponents arising from the Green function we obtain the contribution to the action linear 
in the external field $: 

T dp 



S c =i dt -JL<f,(- p ,t)K c (p,t), (41) 

JO J-oo 27T 

where the "current" 1Z C (which depends on the chosen configuration) is equal to 

TZ C (p, t) = Tit (p) exp (-i|p|»/ti) + TZ f (p) exp {-i\p\v f (T - tj)) , (42) 

and 

^/ (p) = ^2 ( p ) t exp ( ipx ") ~ exp ^P x ')\ + (~P) I ex P (w) ~ ex P (w')] > 

x..;x'..;y..;y>... 

it-i (p)= ^2 ® (~p) t exp ~ exp ^p 5 ')] + (p) t exp (ipy) ~ exp (tpy 1 )] ( 43 ) 

x..;x' ..;y..;y' ... 

for the initial (annihilation operators) and final (creation operators) configurations respectively. Coordinates x, y in 
eq. (|43|l are the coordinates of annihilation and creation operators for the configuration we arc interested in. Finally, 
we get the following functional integral: 

y^^exp I J dtdhj ^( p ,t)^(~p,t 1 )V- 1 (p)S(t-t 1 ) 



1 f°° dt) f^ 1 f°° dv 

-— J dtdhj ^# (-p, i) * (p, ti) |p| exp [-*bl«/l* - *i|] J Q dt J ^®(-P,t)Kc 



(44) 



Here the first term is the action of ea. (|39|l . the second is the quantum determinant and the third comes from the 
Green functions, ea. (|4*TJl . 
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The integral in ea.<]44[l is a Gaussian one: it can be calculated by standard methods. One has to find the saddle 
point field <&o and shift the variables of integration $ — » $ — $ • The integral in fluctuations $ — $ gives the shift 
of the ground state energy due to the electron interaction and normalization coefficient of the ground state wave 
function. We calculate this integral in Appendix [U] The operator structure of the evolution operator is determined 
completely by the terms which appear as a result of substituting the saddle point $o m eq.|@3J. We write them as 
an "effective action" : 



Seff= 2 J ij Q dt ^(p,t)TZ c (p,t). 



(45) 



The saddle point field &o(x, t) obeys the integral equation: 



i T , 1 



* (p,t)-— / dt 1 $ (p,t 1 )\p\exp[-i\p\v f \t-h\]=-iK c (p,t) (46) 
v (p) ^ Jo 

which can be reduced to the following differential equation (to see this it is sufficient to differentiate both sides of 
ea. (|46|l with respect to time): 

d^o(p,t)+u 2 p $ ( P ,t) = 0, (47) 

where 



Wp = |p|«Jl + ^M. (48) 

The boundary conditions for this equation follow from the original integral equation (|46|l : 

a t $ (p, 0) - i\p\v f $ ( P , 0) = 2i\p\v f V (p) Ki (p) 
a t $ (p, T) + i\p\v f $o (p, T) = -2%\p\v f V (p) ft, (p) . (49) 

In the derivation of ea. l|47|l we have used the fact that our system is electrically neutral and hence: 

ft/(p = 0,f) = ft(p = 0,f) = 0. 

The solution of the differential equation for the saddle point field [ea. (|47|l ] gives: 

$o (p, t) = ~ 2 } P } Vf ^ < ; P) c2 , [Ki (exp (-iu p t) + ^ exp (-iuj p (T - £))) 
K + \P\Vf) (1 - £|) 

+ft/ {i P exp (-iw p t) + exp (-iui p (T - t)))} , (50) 

where 



1 _ h + YM 

£p = ; exp \- iuj p T ) ■ 

1 + Ji^YM 



Substituting the saddle point field into the expression for the effective action (|45|) we obtain finally 

o V V(p) 1 ^ 



x [[ft, (-p) ft, (p) + ft, (-p) ft, (p)] F 2 (p) + 2Fi (p) ft/ (-p) ft (p)] , (51) 
where we introduce the following two functions: 

exp (-«|p|u/T) - exp {-iu> p T) 1 - exp (-i (uj p + \p\v f ) T) 

fi (p) n 1- 4p — n 

UJ p - \p\Vf UJp + \p\Vf 

l-cxp(-i(u}p + \p\v f )T) exp(-i\p\v f T)-exp(~iuj p T) 

■ t 2{P) = — j— j h£p j— j • [OZ) 

UJ p + \p\Vf UJp - \p\Vf 
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We return in the expression for the effective action (|51() to a sum over the particle momentum p n 
to the ordinary rule |35j : 



2im/L according 



dp 
2^ 



v 



This will allow us to qualify different infrared divergences which appear in the effective action. Let us note that there 
is no term with p = in these sums. This fact is related to the gauge invariance of the system: constant (in space) 
fields correspond to a pure gauge electrical potential and should not contribute. 

Let us proceed with the wave function of the ground state in Tomonaga-Luttinger model. As was mentioned above, 
in order to separate the ground state we have to take the limit T — > oo. (This corresponds to the case of zero 
temperature.) We can omit oscillating exponents in this limit. As a result we are left only with the function Fi which 
is changed to 




F 2 (p) ~ 

Effective action factorizes into contribution of initial and final states: 

V( P ) 



S, 



eff 



p#o \p\vf 



1 



1 



v(p) 

■KVf 



[Kf(-p)Kf (p) + n i (-p)K l (p)} = 



S eff+ S lff' 



(53) 



Besides S e ff we have to calculate preexponential factors which arise from free Feynman Green functions. At 
T — ► oo only Green functions with equal time arguments survive. As a result, we see that the whole expression for 
the evolution operator for large T factorizes into the product of the wave function of the ground state |fi > and its 
complex conjugate (see also below, section HV Ajl . The final expression for the wave function is of the form 



oo 



n=0 



dxdx' a R (x) b R (x') 



2ni 



iS 



dydy' 4 (y)b[ (y') 
2iri y — y' — iS 



s lff ( x > x> , ~,y,y', ■■) \ F > 



expi>' eff 



(54) 



Let us check, first of all, that the wave function of noninteracting fermions (V = 0) is \F >. The general term in 
the sum of ea. 1)54(1 is a product of factors: 



dxdx' a R (x) b R (x') 



2-rti 



x — iS 



F > 



Let us note now that all singularities of the operator P R [x') are in the upper semiplane (see definition|3Il as well as the 
pole of the Green function. One can close the contour of x' in the lower semiplane and prove that the corresponding 
integral vanishes. The only term which survives is the term with n = and hence |f2 >= \F > as it should be for 
non-interacting fermions. 

A non-trivial answer for the wave function appears only owing to singularities of the effective action. It is clear 
from the general structure of the action (which is the product of Rf(p)Rf(—p)) that the wave function contains only 
terms where both R- and L-particles are present. All terms with only R (or only L) electrons or holes vanish. The 
structure of the ground state wave function is a bit different for different potentials. We will describe it below. 



A. Short range potential. 

We begin with the short range potential: V (p) — Vo- The simplest possible contribution to the ground state wave 
function |fi > (see ea.(|54fl) is 

dxdx' dydy' a R (x) b R (x r ) a{ (y) b[ (y r ) f , 

~^T~- — — / } rj-expS^Ax^^^). (55) 

2m 2m x' - x — id y — y' — id JJ 
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The effective action S e f / for this term has the form: 

S f eff (x, x', y, y') = ^2 — [expip n (x-y + id) + 
Pn >o Pn 

exp ip n (x — y + id) — exp ip n (x — y — \-id) — exp ip n (x — y + id)] , (56) 

where 

V 



v f 



l+y/l + Vo/wVf 



(57) 



The sums in ea. (|ST))l can be easily calculated. We obtain: 



Sl ff (*, y, y>) = ° log jr^tlu"^!^ - < 58 > 
11 tt [x' — y + id) [x ~ y' + id) 

Expression l|55|l describes the simplest possible complex in the vacuum of the interacting fermions. This complex 
has all quantum numbers equal to zero. In fact, it describes electron-electron scattering (in the cross channel). 
Correspondingly, all coordinates x, x' , y, y' are close to each other. In general, this complex does not break down any 
continuous symmetry. 

However in the Tomonaga-Luttinger model a special situation arises. The main contribution to the term i|58|) comes 
from the region x' — y; x — y' — * (of order of the transverse size of the channel) but x — y and x' — y' can be arbitrarily 
large. In other words, the complex decays into RL and _Ri-pairs. As we shall see such a wave function leads to a 
spontaneous breakdown of chiral symmetry. 

Let us consider first the limit of the strong interaction: 



Vn 

— » 1. (59) 
irvf 

In this limit cx/tt — ► 1. It can be seen now that for a/n = 1 the poles (x = x' and y = y') corresponding to free 
fermions are cancelled completely by the fermion-fermion interaction (described by exp(5 e //)) with the effective action 
of ea. H58|) . Instead we obtain new poles in the points x' = y — iS and y' = x + i6. Recalling that b R (a/) is analytical in 
the lower and V L (y 1 ) in the upper semiplane we can integrate further over x' and y' . As a result we get the following 
contribution to the ground state wave function: 



dxa\ (x) b\ (x) J dya\ (y) b\ (y) . (60) 

Thus the 4-particle complex decays into 2 non- interacting " bosons" . They are neutral in the electric charge but have 
a non-zero chirality ±2. 

One can check that no other connected complexes appear in the limit of strong interaction. Let us consider, for 
example, charged complexes. The four fermion contribution is exhausted by ea. (|60|l . so we have to consider a 6-fermion 
complex: 



4 (x) b\ (xQ 4 ( Xl ) tf R (x[) aj (y) bj (y') 
x' — x — id x' x — X\ — iS y — y' — iS 
(x-y + iS) (x!-y + iS) (x' - y' + id) (x[ - y' + iS) ^ 



{x — y' + id) (xi — y' + iS) (x' — y + iS) (x' 1 — y + id) 

This complex, indeed, decays into 2 fermions with x\ — > y' — > x and x' — > y — > x[ (relative distance x — x' supposed 
to be large). These fermions are of the form d^ R (x) d R (x) V L (x) and a^ L (x') b R (x r ) b R (x 1 ). Hence this contribution 
is zero owing to the Pauli principle. One can consider also more complicated configurations which could produce 
charged connected complexes and check that they do not appear in the wave function of the ground state. 
The Pauli principle allows one more complex which describe scattering of chiral pairs: 

a R (x) b[ (x) a\ (x) b R (x) . 
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One can extract the corresponding contribution from the connected part of the general expression eq.j55j). The 
integral over x 1 and y' is easily calculated and we obtain: 

/ dxdyd R (x) b{ (x) d\ (y) P R (y) $(x-y), 

where 

* / \ —iS ( '„ iS 

$(x-y) = 



y — x — 2i5 \ y — x — 2i8 

The function <E> (x — y) is finite at any x, y (even in the point x = y) so its contribution to the integral vanishes in 
the limit 8 — > +0. In other words, in the limit of infinitely strong interaction the chiral pairs do not interact. This 
interaction appears, however, in the next approximations in the inverse coupling constant (see Section llV A lfl . 

To obtain the complete expression for the ground state wave function we have to consider complexes with 
8,12. .. particles and separate the connected parts out of these complexes, this is not necessary, however, since, 
according to a general theorem 36] , the complete wave function is the exponent of the connected complexes |37| and 
we have proved that the only connected complexes are the chiral pairs H6U|I . On the other hand the total chirality C 
of | O > should be zero and only terms with C — can appear in the expansion of |fi >. To take this into account we 
introduce the projector Pc=o onto the state with chirality zero. Then the wave function can be written as 



\n >= VZoPc^oexp 



dxd R (x) b{ (x) + I dyd\ (y) b R (y) 



\F > (62) 



The normalization coefficient Zq is calculated in Appendix[Hl We have already discussed the wave function in Section 
IIIII It corresponds to an unbroken symmetry phase in spite of the presence of an infinite number of chiral pairs with 
zero momentum. If the chiral symmetry is broken, the states with different chirality should be degenerate in energy. 
This is not the case if the size of system is finite — the energy of the state with C = is still minimal and order 
parameter A is zero. 

Wave function Q62JI corresponds to the state with minimal possible energy. Hence it is the wave function of the 
system at 8 = 0. Turning to © ^ let us note that there are two types of exponentials (see below ea. (ISI)l ): first with 
exp(— Vfp/Q) and the second is exp(— vjp/Q) (with renormalized Fermi velocity). Corrections of the second type 
correspond to excitations and we will omit them. However there are no excitations with oj p = Vfp (it can be seen, e.g ., 
by method of bosonization). In fact, these exponentials describe the change of the ground state with temperature [3l|. 
Obviously for O 3> lixvj jL the exponential factor is not small but the preexponential factor, i.e. the Green function 
with imaginary time differences about I/O, gives the smallness. The latter is compensated by the action because it 
is proportional to log I/O. This is the case for the temperature region under consideration. In the opposite case the 
Green function (|31|l is not valid. One can use eq. (|4(J|I but it is impossible to transform the sums over p n to integrals 
in order to get ea. l|31|) . As a result, the Green function will be proportional to a small exponential factor. It can not 
be compensated by the logarithmic divergence from action and the whole term will be small. Therefore the ground 
state wave function (|62|l is valid provided 

O < Qchiral = (63) 

Of course one assumes that the number of states is large , i.e. pfL 3> 1. This allows a transition from the sums to 
integrals in the expressions independent of O (or T). 

In the region of higher temperatures, @ c hirai "C ©, eq.(j2U for the Green function is applicable. In this case after 
same algebraic transformations effective action S e f / (|51|) can be rewritten in the form of 

% = -7En tanh ^et P*/ (-ti n f ^ + ^ t-ti n i (?)] 

^n=£0 lPnl ^ 



where v C f = Vf ^/l + Vo/nvj. If 



2n -r-^ 1 exp „' — exp L 



L n^o lPnl 1 + exp e ' f 



Qchtral < © < 2irv c f /L (65) 
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than equation l|64[l can be transformed to 




~T £ FT exp -^F n f { ~ p) Ui {p) ' (66) 

|39|. So, in the temperature region of interest one should take into account another 4-fermions contribution to the 
ground state: 



dxdx dy'dy' d R (x) a R (x) b\ (y r ) b L (y 1 ) qf 
2ni 2iri x — x + VfT — id y' — y' — VfT + iS 



expS^ ff (x,x,y',y'). (67) 



At lower temperature this contribution is exponentially small. The action for this configuration is 

, (y' - y' - v f T + jS){x-x- v f T + is) 

g (x-y' + i6)(y> -x + iS) ' [ ' 

Thus one has a similar result: a pair a} R (x) b^ L (x) in |f2 > and clr (x) 1>l (x) — in < f2|. However the existence of an 
extra pair implies that the whole chirality C of the state is nonzero. So the states with any C exist. Their energies 
differ by values of the order of 2nvf/L. At the temperature eo. ()65fl these states can be considered as degenerate. Then 
a state with fixed chirality is unstable relative to an infinitesimal interaction which is mixing right and left particles 
(e.g. infinitesimal back scattering). Similar to the theory of superconductivity the real ground state of the system is 
a mixture of states with different chirality, ea. (|24"}> . but a fixed chiral phase 9. 

In order to prove that 

6 C = lu(2tt/L) = 2tvv c } /L (69) 

is the phase transition temperature one should consider the higher temperature region: Q ^> O c . The logarithmic 
contribution to the action S e ff arises from n 3> n m i n ~ LO/27tv c j: 3> 1. At smaller n the logarithmic divergence does 
not exist. So 

s e ff = ~ ri P*/ (-p) n f (?) + n i (-p) n i 0)1 

\n\>n mi „ n 

One can compare it with ea. H64|) . The sums in ea. l|7U|l are calculated in Appendix iDl As a result the logarithms from 
ea. (|S5|l have to be replaced by 

log <^ r^'^j - ), (71 , 

where 



C(G)=^/(9-e c ) (72) 

is the coherence length. 

The r.h.side of ea. (|71|l can be expressed by an exponential integral exponent function with imaginary argument. 
In order to prove that £ (8) is the coherence length one should notice that at the length Ax <g; £ (0) the r.h.s. of 
ea. (|71|l tends to log ((Aa; + iS) / (Ax' + i5)), i.e. in such case the system is characterized by the wave function ea. (|24|l . 
(Indeed, in this case it is possible to repeat the calculations made in Section IjIV A|) if one separates all connected 
complexes by the distances smaller than £ (O).) Thus in a region of a specimen smaller than £ one has a coherent state. 
In the opposite case (distances between pairs, Aa; = \x — y\, are larger than C(©)) expression under the integrand 
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begins to oscillate and the divergence does not exist. As a result, one will have small corrections to the action of 
about 

exp(-*Ax/C(e))C(©) /Ax. 
So the 4-fermion contribution l|67|) leads to a term 

J dxdy 4 (*) &l (x) \F X F\b L (y) a R (y) (73) 

in the evolution operator. Thus at the distance \x — y\ 3> C we have configurations with free bosons. Consequently, at 
such scale the state is noncoherent. Therefore a long range order does not exist at lengths larger than £ (9). One can 
also check this directly. Let us calculate the contribution of the state <|73|) to the order parameter density correlator: 
< \a R (z/i) b\ (j/i) 6^ (x\) aji (xi) \ > . After integration over x\ and t/i one has the contribution of this state to A 2 : 

/>pF)~c><e,. 

As A does not depend on L one has the normal phase (see Section llllf) with low symmetry phase fluctuations. This 
means that C is indeed the phase transition temperature and £ is the coherence length. Besides, one has more 
obvious definition C : 

C(e c )~L. 

In this case the whole system can be described by the broken symmetry wave function Ij24(l . Hence one should think 
that the low symmetry phase is realized if < C . The above discussion should make it clear that this transition is 
smeared over the temperature region about C as it should be for the finite size specimen. 



1. Kosterlitz-Thouless phase. 

Let us prove that a Kosterlitz-Thouless phase 0] is likely to form in the Tomonaga-Luttinger model if one takes 
into account corrections to the action due to irVf/Vo- 

We begin with the case of zero temperature and consider again the 4-fermion contribution eq. (|55|l to the ground 
wave function 

dxdx' dydy' a R (x) b R (x 1 ) a{ (y) b{ (y') \{x-y + i5) (x' - y> + iS) ] a ° > 



2ni 2ni x' — x — iS y — y' — i5 



(x 1 — y + iS) [x — y' + iS) 



where = ce/ir. For simplicity we will consider ao close to unity. Let us consider the configuration with two connected 
chiral complexes separated by a distance R large compared to the transverse size of the channel d: x' — y; x — y' ~ d — > 
and | a; — x'\ ~ R , \y — y'\ ~ R — ► oo. The contribution we are interested in is determined by two cuts: y' = x + i5 
and x' = y + iS and is proportional to 

X dy'x . ,. 1 (x'-yT 



(y'-xf (y'-y) 



The last factor of the integrand is of order the 1/R 1 ao . Distances inside the pair y' — x, x' — y are of order of d. The 
contribution of the distant chiral pairs to the integral f74*J) is 

2(l-a ) 



dxdya R (x)b[(x)a{(y)bUy)(r^^) ° \F > . (75) 



In the temperature region we can consider also contributions of the states with C ^ to the ground state. The 
simplest contribution comes again from ea. l|67|) and has the form 



/ 



2(1-q ) 

dxdy [j^Z^] a R (x)b[(x)\F><F\b L (y)a R (y) (7(i) 
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As can be seen from eq.(J7SJ) and ea. (|76|l the probability to find chiral pairs at the distance R is P(R) = \&(R)\ 2 ~ 
1 / J2 2 ( 1— "o) [3^ . This probability decreases with R but much more slowly than in the theory without interaction. The 
average distance between pairs 

< R >= / dRR P(R) ~ L 2a ° 
Jo 

diverges as L — > 00. It is instructive to consider the same quantities in the theory with non-interacting electrons. Here 

the probability to find a chiral pair is P/ ree = ( \ x -y\ ) ( see ec l-C30- As we have seen it results in the independence 
of A on L. The other limiting case is the system with non-zero density of the order parameter. Here the probability 
to find a chiral pair does not depend on the distance R at all while A oc L. The probability under discussion has 
an intermediate behavior. As a result A will increase with L but the power will be smaller than unity. Both these 
properties can be considered as a definition of Kosterlitz-Thouless phase. 

The temperature of the phase transition O c in the Kosterlitz-Thouless system, at ao < 1, to the unbroken phase is 
of the same order as in the limit of the infinitely strong interaction. Indeed our estimate of C in section TlV Al was 
based on the logarithmic divergence of the action. This divergence exists also at ao < 1 and hence our expressions 
for 9 C and correlation length £ are valid in this case too. 

The wave function for the Kosterlitz-Thouless phase does not have the simple form of eq. I|24l) since the interaction 
of chiral pairs is non-zero. Also chiral complexes with more than two particles are present in the wave function of 
the ground state. However properties of this phase are quite similar to properties of the phase with broken symmetry 
which appears in the limit of infinitely strong interaction. 



B. Cooper channel. 

Let us discuss in brief the case of the attractive short range potential V = VqS(x), Vq < (Gorkov's potential [13 )• 
This potential can be used only if the interaction is sufficiently weak: 

M < 1, (77) 

TTVf 

otherwise the spectrum of excitations acquires an imaginary part (see eq.(0SJ)). This demonstrates that the system of 
one dimensional electrons tends to collapse in this case and the point-like potential has to be modified. This effect is 
similar to the well-known instability in the system of interacting oscillators |40| . 

If inequality l)77|) is fulfilled we can use ea.(|74fl derived above but with ao < 0. In this case the system is in 
the Kosterlitz-Thouless phase too. To prove this statement we consider again four-fermion contribution to the wave 
function with two Cooper pairs separated by distance R (x — y; x' — y' — * 0, \x — x'\ ~ \y — y'\ ~ R ~ L). The 
corresponding contribution to the wave function is of the form: 

al (x)a[ (x)bl (y')b{ (y') , 

and the probability to find such a configuration behaves as i? 2 ( 1_ l a ol), j e decays more slowly than in the theory 
without interaction. 



C. Coulomb interaction. 



We are interested in the long range Coulomb potential also in the limit of strong interaction. However as was 
already mentioned in Section [H] since the length of the three-dimensional screening should be still larger than the size 
of the channel, we must have p/a& 3> 1 [2(|. This means that the parameter related to the strength of the Coulomb 
interaction should satisfy the inequality: 

V(p) 2 2p f 

= log > 1. 78 

TTVf irpfcn, p 

This can be large only due to a large logarithm. The argument of this logarithm is of order pfL/n 3> 1. For this 
reason the integrals over momentum in the expression for the action S e ff, ea. (|53|l . should be cut at 

/— irpfcit 1 \ 

Pmax ~ tyf exp I I <C p f 
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This cut-off is not essential if p ma xL S> 1. Then every term of the action can be written as 

— exp (ipSx)/ (l + 3/2 v / 7 rw / /e 2 log~ 1/2 (2 Pf /p) 

l/L P V V 



log Sx — 6y-KVf/e 2 log 1 / 2 (2pf5x) . 

Separating the four-fermion contribution in two chiral pairs a R (x) b L (x) a L (y) b R (y) , we find that the probability to 
find the chiral pair at a large distance \x — y\ from its counterpart behaves like 



exp 



-12yjnv f /e 2 log 1 12 (2p f \x - y\) 



It decays more slowly than any power of \x — y\ 41] . Strictly speaking, this behavior does not correspond to the 
Kosterlitz-Thouless phase but it is clear that physically these two phases are quite similar. Let us note also that in 
the formal limit (nvf/e 2 ) log(2p/L) — > we have again the condensate of independent chiral pairs with the wave 
function of ea. (|24f> . 
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APPENDIX A: EVOLUTION OPERATOR FOR SYSTEMS OF FERMIONS. 

In this appendix we derive the representation for the evolution operator of fermions in the external field as a 
functional integral with definite boundary conditions. 

In the Schrodinger representation the evolution operator S(T) is 



S[T] = Texp(-i [ Hdt)\F >< F\ 
Jo 



where H is the fcrmion Hamiltonian in the external field which is bilinear in the fermion fields. As we have seen the 
general problem with electron- electron interaction can be reduced to this problem at the price of the integration over 
the external field. For simplicity we begin with a model with an empty ground state |0 > rather than the Fermi one. 
(This allows us to write the equations in a more compact form.) Besides we will omit the spatial arguments. 

Let us divide the time interval T in N infinitesimal pieces 5 — T/N (with the point i = N corresponding to t = 
and i = 1 to t = T) and introduce the sum over the complete set of quantum mechanical states \k >< k\ in the 
intermediate points: 

S[T] = \k N >< k N \ (1 - iSH) \k N ^ >...<k 2 \(l- iSB) \ki >< h\. (Al) 

n j 

For any complete set of wave functions in the secondary quantization representation we have: 

exp (Tr&a f ) |0 >< 0| exp (-Tr^a) , (A2) 
(index n corresponds to the set of all quantum numbers). The Grassmann variables £ are defined in the usual way: 
J d^{n) =0,J dti(n)ti(n) = 1, [£(n),&(n)] + = 0.,P^ = J[d^(n). 



Ea. (|A2p can be proved by direct comparison of the left and right-hand sides. We use this representation to rewrite 
the sum over states as a functional integral. 
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At every point i we obtain the following matrix element of the Hamiltonian: 

exp (-Tr^fi) < 0| exp (-Tr^a) (l - 18H (a\a)) exp (-Tr^+^+i) (A3) 

exp (Tre j+ ia f ) |0 > . 

To calculate this matrix element we move all creation operators to the right. For the Hamiltonian H depending 
linearly on a, a + e.g. for the Hamiltonian in the external field the result is 

exp (Vr£l - &) + iSTrH (ej, 6+1)) • 

Thus the result of the calculation is that creation-annihilation operators in the Hamiltonian are substituted by 
Grassmann variables 

The product over all intermediate points at N — > 00 tends to 

exp(- [ dt ¥ (t) [St + iH] * (t)) = expi [ dtjC, C = * [id t -H]® 
Jo Jo 

where £ is the Lagrangian of the system. This expression should be integrated over variables ^ at all intermediate 
points in time. The boundary points are specific, however. Creation operators entering \k^ > and annihilation 
operators entering < k± | are not contracted. They are variables on which the evolution operator depends. 

Let us integrate over all intermediate variables and consider the answer as a function of Grassmann variable 
(and £ N ). This function can be only linear one: A + B^i. Then the last integration in is of the form: 

J V&Vtl exp (-Tr£}&) exp (-Tr^a) (A 1 + TrBifc) = A 1 + TrB x a 

Thus, we see that the variable £1 should be substituted by an annihilation operator. Integrating over £jy we conclude 
that £jy is substituted by a creation operator. 

Finally, we can formulate the following recipe: to calculate the evolution operator one has to integrate exp(i C) 
over 'J', ^ at all intermediate points. At t = the value of VP is fixed to a, at t = T ^> is fixed to a + . The values of ^> 
at t = and ^ at t = T remain arbitrary. As a result, operators a and d + arc defined at different times. Therefore 
one should consider theirs here as anticommutating. 

If the ground state of our system is a filled Fermi sphere we have to introduce two type of creation-annihilation 
operators 0/^,6^ corresponding to electrons and holes. Then we can apply the above derivation to this case as well. 
One should introduce negative (^~) and positive ( V I' + ) frequency parts of W variables and double the number of 
variables £. 



APPENDIX B: CALCULATION OF Det$. 



We will calculate the functional integral over the fields \ an d X- They obey zero initial conditions. 



Det<S> = J VxV X exp(^J^ dt J dx% {idt - n ext (x)) 



(Bl) 



where H ex t — Ho (x) + <& (x, t) . In an ordinary case Det<& can be calculated in the usual way using the identity 

log [Det<S>] = Trlog$. 
After differentiation over A the right-hand side of this identity is represented as 



if d\(-d t -iH (x) -i\®(x)) 1 $ (x) 



Sp 



where the inverse operator is the Green function with the same arguments. Usually the result does not depend on the 
order of arguments. However in the theory with Adler-Schwinger anomaly the sequence of time and spatial arguments 
is essential. The simplest way is to assume to make spatial arguments equal first. In this case the result will be in 
contradiction with gauge invariance of the theory. In the paper 20] the procedure free of this difficulty was suggested. 
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It does not exist in this case because all calculations are done with non-equal variables up to the end. The procedure 
is based on the Heisenberg equation for the electron evolution operator S ($) = Det& exp Sq (T) at the external field 
(without direct electron-electron interaction) . In Heisenberg representation one has 



.dS_ 



H ex t, S 



where H ext is the noninteracting electrons Hamiltonian (the external field is dependent on the time T), and the 
action Sq is defined by the ea. (|34|l . One should note that all creation operators are defined at the moment T and the 
annihilation operators at t = 0, therefore in the commutator [H ext ,So (T)] only the terms with creation operators do 
not commute with Sq. 

One can rewrite the latter equation in the following form: 



.<91ogDei$ 
' df 



exp (—<S ) [Hext, exp So] - 



(B2) 



In order to calculate the commutator in this equation one can use the well-known identity: 



a (x) , exp dx'K (a/) (fi (x r 
= J d Xl A (xi - x') K (xi) exp (J dx'K (x') d f (x') 



which can be proved by expanding the exponential functions. (Here K is anticommutating with d operator and 
A (xi — x') - is the anticommutator {d (x) , at (x')} is defined by eg. 110(1 '). The left-hand side of ea. (|B2l) is a c-number; 
this means that all operators from right-hand side of this equation have to vanish. The c-number parts arise only 
from the following commutators: 



dx$ (x) 



b (x) d (x) , exp / dydy'tf (y') G (y'T, yT - e) S+ (y) 



As a result one has 



.d log Det$ 



dxdydy' 

2-$(x,T) 



[W 



dT J (2ni) 

G L {y'T lV T-e) 

iS) (y — x — iS) (x — y' — iS) (x — y — iS) 



G R (y'T,yT- 



(B3) 



This representation is general. In order to rewrite the right-hand side of this equation in our case one should take into 
account that only the region y — > y' — > x is essential in the first term. However at the point y — > y' the argument of the 
exponential in the Green function (|32|l vanishes. This means that the contribution is determined by the preexponential 
pole and only the first and the second terms of the expansion of exponential can give nonvanishing contributions. All 
singularities under the integral over y in the function coming from the first term are in one semiplane. One can close 
the contour in the other one and prove that this integral vanishes. In the next order in <f> only the part with the 
singularity in the lower semiplane of y gives a nonvanishing term. After integration over y' one has (in momentum 
space representation ): 



-i 
2n 



7T- / dtl 



^P$- P (T) $ p (ii) exp (-ipv f (T - h)). 
Zn 



The L-electrons give the same result but with the opposite sign of p in the region p < 0. After the integration of 
ea. (IB3(l and symmetrization one gets ea. (|35ll . Note that ea. 1)35(1 is gauge invariant: the fields depending only on time 
do not contribute to ea. ((35(l . 



APPENDIX C: NORMALIZATION COEFFICIENT AND ENERGY SHIFT. 



We have seen that the matrix element from ea. l(13|) can be expressed as a Gaussian type functional integral. It 
gives the normalization coefficient and the ground state energy shift. Indeed, one can expand the exact wave function 
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over the free electrons ones. In the limit T — * oo only the matrix element between the lowest energy level survive. It 
can be represented in the form: 

Z = exp(~iAET)\ < fl\F > | 2 

where AE is the ground state energy shift. Comparing Z with the definition of the normalization coefficient Zq from 
eq. (|24J) one can see that it is equal to the overlap probability of the ground states of the free and interacting electrons : 
| < fl\F > | 2 . The normalization coefficient should be calculated for a finite size system, as it is exponentially small 
with the volume. 

On the other hand, the matrix element we are interested in is 



dtdti 



^(p,t)^(- P ,t 1 )v- 1 (p)Ht-ti) 



Z = 1/N I P$exp 

1 f 1 f°° dv 

dtdh /$(- P ,i)$(p,OI;p|expHI;PM*-<i 



(CI) 



where 1/7V is the normalization coefficient l|38JI ■ 

It is convenient to transform integral operator The more effective procedure is the transformation of the integral 
operator IClfl to a differential form. In order to do this we note the identity 



-2i 



Thus symbolically 



exp(-i\p\v f \t-ti\) 
and the kernel of ea. (|Cl|) is equal to: 



1 ( d 2 \ f T 

Wf \W +pV U J ^( f i) ex PHH"/l i - i il) = ,1, p( i ) 

-2i\p\vf 



(C2) 



9 4- , , 2 



■ p 2 vj 



6(t-h). 



As a result, one has: 

Z = l/JV I 2?$exp 



dt 



Z7T — — -L -n z ')'^ 



or 2 



p*v 



(C3) 



(C4) 



Taking into account that the normalization coefficient Af cancels the sjDet (i/2V p ) that arises from the differential 
kernel definition one gets: 



Z- A = Det 



D 

Do" 



In order to define the differential operator one should have two initial conditions. In exactly the same way as for the 
derivation of the equation for the saddle point field one gets this conditions: 



dt$(p,0)-i\p\vf<S>(p,0) = 
d t $(p,T)+i\p\v f $(j3,T) = 0. 

Usually, one calculates determinants with zero boundary conditions 

% (o) = $ p (T) = o 

In order to reduce our problem to the problem with zero boundary conditions let us introduce 

% (t) = ¥ p (t) + <$> {t) . 



(C5) 
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The field $ p (t) is supposed to obey the equation A t $ p (t) = 0, with initial conditions ea. (|C5(l . (As usual A t = 
^2 + <-^p-) The field (t) is arbitrary but with zero boundary conditions. The solution is $ p (t) can be expressed in 



the following form: 



— / \ , , sinujr, (T — t) sinw„(t) 

<E p (t) = $ p (0) / V - > +<t> p (T)- ■ 



sin w p T 



sin uipT 



(Constants <J> P (0) and $ p (T) are arbitrary.) This means that the determinant has the form 



D 



-1/2 _ 



d% (0) d% (T) exp [j ($ p (T) a t p (T) - $ p (0) (0))] 

V<f> p (t) exp (^_ p (t) A^p(t)). 



The integral over ^ p can be calculated in the usual way: 

C(p = 0) 



sin w p T ' 



(C6) 



where C (p = 0) is the p = contribution. It will cancel in the final expression. Taking into account the identity 
$ p (T) d t ((> P (T) - $ p (0) 9 t p (0) = - ($ p (0) 2 + $ p (T) 2 ) (|p|t;/ - iw p cot u p T) + 



sin u>„T 



^-$ p (0)$ p (T) , 



one has 



z -=n 



sin 



LOpTp vj — 2iLu p \p\vf cotui p T + uj. 



p^O 



sin|p|T 2\p\vfLu p (l — icotujpT) 



(C7) 



If temperature is non-zero one should substitute T by 1/0. Note that ea. l|C7|) is valid even at the temperature 
<C QcMraU because only the Green functions with equal time arguments were used. In this temperature region the 
Z can be expressed in the following form: 



z = n exp 

pjtO 



26 / uj p + \p\vf 



(C8) 



It is convenient to rewrite this equation in the following way: 



Z = exp 



Jo 2tt ^ w p 



4-\/| pk/^p 



(C9) 



In this form one can see the energy shift (first term in the exponent)and normalization coefficient (second term) 
explicitly. The sums in this equation diverge because of the gapless spectrum. They have to cut off at p ma x ~ 1 1 d. In 
order to take the preexponential factor one should calculate the next correction after the Riemann sum. As a result 
we have for the short range potential 



for the energy shift and 



Z a = 4W-TO//Vbexp 



And 



\ogV /irvf 



for the normalization coefficient. 
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APPENDIX D: CALCULATION OF SUMS IN SECTION IIV Al 



All the sums in the equation for the action can be calculated by differentiation of the expression S (a) with respect 
to parameter a: 



5(a) 



2tt ^ 1 

-p }. — exp 



2mna 



(x + iS) 



( a varies within the region (1, zoo)). After summation of the geometric series one can rewrite it in the following form: 

- l - ya dy 



S(l,x) 



y 



(i-yY 



where yo (x) = exp (2ni/ L) (x + iS) . It is implied that x <C L here. This result can apply at x ~ L as an order of 
magnitude estimate only. The final expression appearing in the action is: 



S(l,x)-S(l,y) = 



{x+lS) dz 
(v+iS) z 



exp 



iz 

T 



(Dl) 



where £ = L/2tt (n min — 1)~ .If n min ~ LQ/2nvj then ( is equal to the coherence length ea. (|72|) . 

Let us consider the influence of the boundary conditions on the action. In principle, any of them can be rewritten 
as p n = 2ir (n + Sn) /L, \5n\ < 1/2. In this case at 6 = the action is determined by the sum: 



S' (a) = exp 



2ni (n + Sn) a 



(x + id) 



In the same way one gets 



S'(l,x)-S'(l,y) 



{x+ ' s) dz 
(y+iS) z 



■ exp 



2-Kizbn 
L 



The result of this is that up to \x — y\ ~ L at = the action does not depend on the boundary conditions. For 
(0 3> C ) one should cut off the sum at some n = n m i n . As a result £ is substituted by L/2n (n m i n + Sn — 1) in 
ea. (|Dl|l . This suggests the replacement of C by (1 — Sn) C 01 • However the transition temperature can be defined 
only up to a factor of the order of unity. Therefore we should not take this into account. 
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